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ABSTRACT 

We present results from Chandra and XMM— Newton observations of the radio pulsar B1046— 58. 
A high-resolution spatial analysis reveals an asymmetric pulsar wind nebula (PWN) ~6"xll" in size. 
The combined emission from the pulsar and its PWN is faint, with a best-fit power-law photon index 
of r=1.7^g'2 an d unabsorbed luminosity of ~10 32 ergs s _1 in the 0.5—10.0 keV range (assuming a 
distance of 2.7 kpc). A spatially resolved imaging analysis suggests the presence of softer emission from 
the pulsar. No pulsations are detected from PSR B1046— 58; assuming a worst-case sinusoidal pulse 
profile, we derive a 3tr upper limit for the pulsed fraction in the 0.5-10.0 keV range of 53%. Extended 
PWN emission is seen within 2" of the pulsar; the additional structures are highly asymmetric and 
extend predominantly to the south-east. We discuss the emission from the PWN as resulting from 
material downstream of the wind termination shock, as outflow from the pulsar or as structures 
confined by a high space velocity. The first two interpretations imply equipartition fields in the 
observed structures of > 40—100 /j,G, while the latter case implies a velocity for the pulsar of > 190 

Uq 1 ^ 2 km s~ 4 (where no is the ambient number density in units of cm -3 ). No emission from an 
associated supernova remnant is detected. 

Subject headings: pulsars: general — pulsars: individual (PSR B1046— 58) — X-rays: general 



1. INTRODUCTION 

X-ray observations of rotation-powered pulsars repre- 
sent a powerful tool for studying the energetics and emis- 
sion mechanisms of these objects. A large fraction of the 
available rotational energy is thought to be carried away 
in a relativistic wind of particles. When this wind is 
confined by the surrounding medium it decelerates and 
a synchrotron nebula forms, called a pulsar wind neb- 
ula (PWN) . The overall PWN characteristics provide in- 
sights into the content and energy spectrum of the pulsar 
wind, the large-scale magnetic fields, and the surround- 
ing medium. The small-scale structures of the nebula 
provide details of the acceleration sites, instabilities in 
the magnetic field, and the interaction between the wind 
and its surroundings. In young, energetic pulsars, non- 
thermal emission from particles accelerated in the mag- 
netosphere and thermal emission from the surface of the 
star are also expected to be present. 

The radio pulsar (PSR) B1046— 58 was discovered dur- 
ing a Parkes survey of the Galactic plane ({Johnston et al.l 
1992). It has a period of P = 124 ms and period deriva- 
tive P — 9.6xl0 -14 s s _1 . These values imply a char- 
acteristic age of t c = P/2P = 20.4 kyr, spin-down lu- 
minosity of E = 4ir 2 IP/P 3 = 2.0xl0 36 ergs s" 1 (with 
a moment of inertia of the neutron star / = 10 45 g 
cm~ 2 ) and surface dipole magnetic field strength of B = 

3.2xl0 19 VPP G = 3.5xl0 12 G. These properties are 
similar to those of other young neutron stars typified 
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by the Vela pulsar. The dispersion measure of 129 pc 
cm" 3 towards the pul sar implies a distance of 2.7 kpc 
(jCordes k Lazioll2002D . 

Deep radio observations did not detect extended emis- 
sion associated with PSR B1046— 58 ijStappers et alJ 
1999). X-ray observations with ASCA and ROSAT de- 
tected emission near the pulsar and suggested the pres- 
ence of large-scale structures surrounded by faint emis- 
sion ijPivovaroff et aUl2000a|) . However, the poor angu- 
lar resolution of these observations prevented conclusive 
interpretation of the data. PSR B1046— 58 is also one 
of a few pulsars with a possible EGRET 7-ray counter- 
part. The EGRET source 2EG J1049-5827is coincident 
with the radio coordinates of PSR B1046— 58 and likely 
7-ray pulsations a t the radio period have been found 
ffCaspi et al.ll200rH) . 

Here we report on observations of PSR B1046— 58 car- 
ried out with the Chandra and XMM— Newton satel- 
lites. Our study reveals a faint, arcsecond-scale PWN 
surrounding the pulsar. A detailed imaging and spec- 
tral analysis of the system suggests the presence of soft 
emission from the pulsar. We also examine the charac- 
teristics of the nebula in light of current models for the 
production of PWNe. 

2. OBSERVATION AND DATA ANALYSIS 

PSR B1046-58 was observed with XMM-Newton on 
2002 August 10. The European Phot on Imaging Camera 
(EPIC) MPS and P N instruments ijTurner et alJl200H 
IStriider et alJ I2001|) were operated in full- and small- 
frame modes, respectively. These settings provide a tem- 
poral resolution of 2.6 s for EPIC-MOS and 6 ms for 
EPIC-PN. The medium filters were used for MOS, the 
thin filter for PN. The CCD data were reduced with the 
XMM-Newton Science Analysis System (SAS v6.2.0). 
The standard screening criteria results in an exposure 
time for the MOS cameras of 20 ks. For the PN camera, 
only ^70% of the exposure time was used in the small- 
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window mode (4' field of view), resulting in an effective 
exposure time of 15 ks. 

The Chandra observation of PSR B1046— 58 was car- 
ried out on 2003 August 8 -9. The aimpoint of t he back- 
illuminated ACIS-S3 chip ijGarmire et al]l2003D was po- 
sitioned at the radio coordinates of the pulsar. The ob- 
servation was taken in timed exposure, very faint mode, 
providing a temporal resolution of 3.2 s. The data were 
reduced using the CIAO 3.2.0 software and standard rou- 
tines. The resulting effective exposure time was 36 ks. 

The EPIC-PN instrument offers a higher sensitivity 
than ACIS-S, but its lower exposure time and reduced 
field of view limited its use to spectral and temporal 
analyses of the pulsar emission. In turn, the Chandra 
data were used to perform spectral and high-resolution 
imaging analyses of PSR B1046— 58 and its surroundings, 
while the EPIC-MOS observations were mainly used as a 
consistency check due to their lower sensitivity and lower 
spatial resolution. 

3. IMAGING ANALYSIS 

Figure ^ shows the area surrounding the radio posi- 
tion of PSR B1046-58 obtained with Chandra. The im- 
age has been smoothed with a Gaussian with er=0'.'5. 
It reveals for the first time the detailed extended struc- 
tures surrounding PSR B1046— 58, which we designate 
as its pulsar wind nebula (PWN) based on its over- 
all characteristics discussed in the next sections. The 
nebula is elongated and ~6"xll" in size, with its ma- 
jor axis oriented in a SE-NW direction. The line in 
Figure ^ marks the radio coordinates of the pulsar at 
a^ nn ^10 /l 48 m 12!6(4"7) and ^n nn=-58 o 32 / 03 / /75(0 , /02), 
2cr errors (Sta ppers etaD 11999). While there is emis- 
sion immediately surrounding the pulsar position (PWN 
"head" ) , there is also clumped emission to the SE (PWN 
"body" ) and a bright clump to the NW ( "north clump" ) . 

The head region of the nebula is consistent with a 
point source surrounded by extended structures. Fig- 
ure |2 shows the intensity profile of the observed emission 
along a 10"x35" region in a SE-NW direction aligned 
with the implied axis of symmetry of the nebula. The 
shaded region represents the 2a background level ob- 
tained from a nearby, source-free region. The dashed 
line represents Chandra's PSF at an energy of 1.5 keV 
obtained using CALDBv3.0.3 and the mkpsf tool 7 . As 
the head region is also consistent with the radio coor- 
dinates of the pulsar, we consider the source embedded 
in the region to be the X-ray counterpart of the pul- 
sar. Comparing the X-ray position of other sources in 
the field with their optical counterparts we derive an X- 
ray position for the pulsar of a2ooo=10 ?i 48 m 12?64 and 
<*>20oo= — 58°32'03'/6, with an overall rms error of 0'.'55. 
This is coincident within la of the radio position, espe- 
cially in the more constraining declination direction (see 
Figure We then define the "pulsar" emission as that 
arising in a circle of 1" in radius centered on the above 
X-ray coordinates. 



7 The slight asymmetry in the PSF profile is due to the bin- 
ning and extraction region used in order to match the pro- 
file settings for the pulsar. A small, intrinsic asymmetry in 
the PSF has been suggested and could also contribute (see 
|http: //cxc. harvard. edu/cal/docs/cal_present_status.h tml#psf I. 
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Fig. 1. — Combined Chandra image of PSR B1046-58 and 
its PWN using individual images in the 0.5—2.0 keV (red) and 
2.0—10.0 keV (blue) ranges. Each image was smoothed with a 
Gaussian with <r=0"5. The circle (0"55-radius) and line (4'.'7 
length) mark the X-ray and radio positions of the pulsar, respec- 
tively (see §3). The main components discussed in the text are 
labeled. 
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Fig. 2.— Chandra axial profile of PSR B1046-58 and its PWN 
in the 0.5—10.0 keV range in a SE-NW direction (solid line) aligned 
with the implied axis of symmetry (see Fig. 0. The origin repre- 
sents the X-ray position of the pulsar. The extended emission at 
offsets of ~— 3" to —8" arises from the "body" of the PWN (in the 
SE direction) while the emission at offsets of ~3"— 5" arises from 
the "north clump". We also show the 2a level of the background 
(shaded region) along with the PSF of the telescope at 1.5 keV 
(dashed line). 
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Fig. 3. — Chandra image of the field surrounding PSR B1046— 58 
obtained from the ACIS-S3 chip and part of ACIS-S2. Individ- 
ual background-subtracted and exposure corrected images in the 
0.5-2.0 keV (red) and 2.0-10.0 keV (blue) ranges were used. The 
images were smoothed with a gaussian with cr=2" . The circles 
mark the position of the so urces detected with AS CA and their 
uncertainties (25" radius, see IPivovaroff et al. 2000a), all of which 
are coincident with resolved point sources in the present observa- 
tions. Similar results (although at lower angular resolution) were 
obtained with the XMM— Newton MOS instruments. No extended 
PWN or supernova remnant emission was detected in the data. 



We also find that the extended emission from the 
PWN dominates at high X-ray energies. Tabled shows 
the number of background-subtracted counts obtained 
for various parts of the system in the soft (S) and 
hard (H) bands, their hardness ratios (HR = S— H / 
S+H) and their total signal-to-noise ratio (S/N, e.g., 
IPivovaroff et all l2000b(l in the 0.5-10.0 keV range. 
While emission from the nebula is likely to be present 
within the 1" radius used for the pulsar, we suggest that 
emission from the neutron star is indeed present (and 
may dominate) in this region. The above is supported 
by the fact that this region is the only one that is seen to 
exhibit a positive HR with high significance, suggesting 
a contribution of distinctly softer emission than that of 
the rest of the n ebula. 

Using ASCA, IPivovaroff et all ( 2000a) detected four 
point sources within ~4' of the radio position of the 
pulsar embedded in possible faint emission. The Chan- 
dra and XMM— Newton observations confirm the pres- 
ence of these X-ray point sources in the field, in addi- 
tion to many others resolved by Chandra (see Figure |3J). 
The soft point source in Chandr a coincident with "Src 
2" from IPivovaroff et alJ l)2000al) h as an optical coun - 
terpart in the USNO B1.0 catalog ijMonet et al]l2002j) . 
while we find no counterparts for the hard sources co- 
incident with "Src 3" and "Src 4" in present catalogues 
(e.g., USNO B1.0; 2MASS, Cutri et al. 2003; RASS, Vo- 
ges et al. 2001). No additional extended emission was 
detected that can be attributed to the PWN beyond the 
above arc-second structures. We then conclude that the 
low angular resolution and broad PSF wings of the tele- 
scope significantly affected the ASCA observations due 
to the large number of point sources in the field. We also 
note that the ASCA positions appear to be systemati- 
cally north of the Chandra coordinates, the latter being 



Fig. 4. — Chandra spectrum of the combined emission from 
PSR B1046-58 and its PWN. The data were binned to contain a 
minimum of 20 counts bin -1 . The solid line represents the best-fit 
power-law model shown in Table l2*| and the lower panel shows the 
fit residuals. 



in agreement with source positions at other wavelengths. 

4. SPECTRAL ANALYSIS 

The data obtained from both Chandra and 
XMM— Newton were used to perform a spectral analysis 
of the emission from PSR B1046-58 and its PWN. For 
the Chandra data, an elliptical region of 12"xl8" in size 
was used with a surrounding annulus as background. 
A total of 184±14 background-subtracted counts were 
detected in the 0.5-10.0 keV range. The XMM- Newton 
data were extracted from circular regions of 25" radii in 
each detector, encompassing ^77% of the photons from 
a point source, with a nearby region used as background. 
The data from MOS1 and MOS2 were not included in 
the analysis due to the low number of counts detected 
in them (106±17 and 47±17, respectively 8 ). The 
PN detector collected 260±47 background-subtracted 
counts in the same range. The extracted spectra from 
Chandra's S3 chip and XMM- Newton's PN detector 
were fit simultaneously in the 0.5—10.0 keV range using 
XSPEC (v. 11. 3.0) and a minimum of 20 counts per bin 
in each spectrum. 

Thermal models for the integrated emission from 
the PWN and PSR B1046-58 are statistically ac- 
ceptable (xl^lA). However, they result in temper- 
atures far too high to represent shock-heated ther- 
mal plasma (T>8xl0 7 K, Raymond-Smith model, 
IRavmond fc Smitbl Il977l) or emission from the surface 
of the neutron star (T>lxl0 7 K, blackbody model). 
Instead, a non-thermal absorbed power-law model pro- 
duces an equally acceptable fit (see Table |2J| with a pho- 
ton index of F=1.71q'2, similar to what is observed for 
other pulsars and their PWNe (e.g., iGaensler fc Slanel 
2006). The Chandra spectrum with the best-fit power- 
law model is shown in Figure^ The observed emission is 
thus consistent with having a predominantly non-thermal 
origin. 

8 It has been found that the MOS1 count rate can be arti- 
ficially boosted by statistical fluctuations by up to ~70%, es- 
pecially for sources with <50 true counts . See §2.1 and §4 i n 
http: / / xmm.vilspa.esa.es / docs / documents / CAL-TN-0023-2-1 ,ps 
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We then searched for differences in the spectral char- 
acteristics of the different components of the system, as 
is hinted from our imaging analysis. The high spatial 
resolution available with Chandra was required in this 
case. The 116±14 counts detected from the PWN alone 
were fit with a power-law model. Holding the interstellar 
absorption fixed to the best-fit value from Tabled the re- 
sulting photon index was Tpwn ~ 1.2. Holding this con- 
tribution from the nebula fixed while fitting the overall 
spectrum, the residual (presumably pulsar-dominated) 
emission was well described by a power-law model with 
a photon index of T ~ 2.4. If a blackbody model is fitted 
to this emission the resulting temperature is ~6.1xl0 6 
K, too high to represent purely thermal emission from 
the neutron star but lower than thermal fits to the over- 
all spectrum. Although the small number of counts did 
not allow us to constrain the above parameters, emission 
from the pulsar softer than that of the nebula seems to 
be present. 

5. TIMING ANALYSIS 

The data from the XMM- Newton EPIC-PN camera 
were used to search for pulsations from PSR B1046— 58. 
No evidence for an i nstrumental 1-s jum p in the data 
was found (see, e.g.. iWoods et alJ 12004). The dataset 
was converted to the solar system barycenter and a cir- 
cular region of 25" radius centered at the above Chandra 
coordinates was used to extract the pulsar events. We 
also examined the dataset in different energy ranges us- 
ing 0.5-10.0 keV, 0.5-2.0 keV and 2.0-10.0 keV. Radio 
observations of the pulsar obtained with the Parkes tele- 
scope predict a period for the middle of our observation 
(MJD 52496.5) of / = 8.08512306 Hz (P = 123.683955 
ms). 

The Z% test l|Buccheri et al.l fi983) was used to search 
for a periodic signal by folding the extracted photons over 
a range of 10 trial frequencies centered on the radio pre- 
diction. The number of harmonics used was varied to be 
n = 1, 2, 4, and 8. The most significant signal was found 
in the 0.5-2.0 keV range at /=8.08517(4) Hz (la error) 
with Zg = 8.9. Although this signal is consistent with 
the radio period of PSR B1046— 58, it has a probability 
of chance occurrence of 6.3% for a single trial and it is 
not signi ficant given the numbe r of searc hes performed. 
Follow ing I Vaughan et all ljl994j) and, e.g,. IRansom et all 
( 2002), we can derive an upper limit on the pulsed frac- 
tion. Assuming a worse-case sinusoidal modulation, the 
maximum Fourier power obtained for a small range of 
frequencies centered on the radio prediction can be used 
to calculate an upper limit for the pulsed fraction at a 
specific level of confidence. In this way, we derive an up- 
per limit for the pulsed fraction at the 99% confidence 
level of 53%, 65% and 61% in the 0.5-10 keV, 0.5-2.0 
keV and 2.0—10.0 keV ranges, respectively. 

6. DISCUSSION 

6.1. PSR B 1046- 58 

X-ray emission from young radio pulsars is expected to 
include contributions from any of the following processes: 
thermal emission from the entire surface due to initial 
cooling, non-thermal emission from magnetospheric pro- 
cesses, and possibly thermal emission from localized hot 
spots reh eated by back-flow ing magnetospheric currents 
(see, e.g.. lKaspi et al.ll2004|) . 



The combined emission from the pulsar and its PWN is 
best described by a non-thermal power-law model. How- 
ever, our imaging and spectral analysis suggests that an 
additional, soft component may be present in the pulsar's 
emission. The soft emission can be well described by a 
steep power law with T ~ 2.4, although such an interpre- 
tation would contradict the generally observed trend of 
young pulsars havin g spectral indic es harder than those 
of their PWNe fe.g.. lGottheTll2003|) . The additional soft 
emission cannot be entirely thermal, as the derived tem- 
perature is too high. It could, however, represent a com- 
bined spectrum of (hard) non-thermal plus (soft) ther- 
mal emission. The number of counts detected does not 
allow us to fit two-component models to this emission 
alone. Instead, in the case that this overall soft excess is 
related to thermal emission from the whole surface, we 
can constrain its temperature by fitting a blackbody plus 
power-law model and finding the maximum temperature 
that yields an acceptable fit to the total spectrum. We 
adopt hard limits for the interstellar absorption of Nh < 
2x 10 22 cm -2 and a blackbody radius observed at infinity 
of > 10 km (at a distance of 2.7 kpc). The power- 
law index and normalization were allowed to vary freely 
and the maximum blackbody temperature that causes 
deviations from the best-fit model at the 3a level was 
found. This results in a surface temperature of Tg? < 
1.4x 10 6 K. This temperature is not constraining for cool- 
ing models of neutron stars, as the temperature in the 
case of minimal cooling is predicted t o be < 1-1 xlO 6 K 
for the pulsar's characteristic age (e.g-. IPage et all 20051. 
Similarly, if emission was present from a small polar cap 
on the surface with R°° < 3 km, the resulting temper- 
ature is T§g > 2.5xl0 6 K . This is comparab le to those 
seen in other pulsars fe.g., IPavlov et afll2002[) and could 
then explain the observed soft excess for PSR B1046— 58. 

6.2. The PWN 

The high spatial resolution available with Chandra al- 
lowed us to discover the arc-second scale PWN structures 
surrounding the pulsar. The "head" of the nebula is coin- 
cident with the radio position of the pulsar, while diffuse 
emission is seen predominantly to the SE. The emission 
from the pulsar and PWN is very faint, with a com- 
bined unabsorbed luminosity of only ^lxlO 32 ergs s _1 
in the 0.5—10.0 keV range. The efficiency with which the 
pulsar converts its rotational kinetic energy into X-rays 
is then r\x = Lx/E ~ 5xl0~ 5 , comparab le to the val- 
ues found for other Vela-l ike systems (e.g., iPavlov et aD 
12001 ttlCarn"ilo et, a.Ul200l . 

It has been recently suggested that an empirical rela- 
tionship exists between the X-ray spectral power-law in- 
dices of young pulsars (T psr ) and their PWN (T pwn ) with 

the pulsar's spin-down energy loss E l|GotthelJ I2003) . 
According to this relationship, lower E pulsars exhibit 
harder spectral indeces. For PSR B1046— 58, the pre- 
dicted values are r psr ~ 0.1 and T pwn ~ 0.9. The PWN 
emission appears to exhibit a spectrum similar to that 
predicted by this relationship (T pwn ~ 1.2). However, 
the small number of counts and the possible detection of 
thermal emission limit our ability to constrain the spec- 
tral characteristics of the pulsar and its PWN. Obser- 
vations of other Vela-like pulsars do not seem to sup- 
port this relationship (e.g., PSR B1823— 13, Gaensler et 
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al. 2003). However, longer observations at high angular 
resolution are needed to reliably test the validity of the 
relationship in this group of pulsars 9 . 

Radio observations did not detect emission from the 
PWN and constraints on the radio properties depend 
directly on the underlying assumption for the radio ef- 
ficiency (tjr = Lr/E; iStappers et alJll999|) . Detected 
nebulae show efficiencies in the range vr ~ 10~ 4 — 10~ 3 , 
while up per limits for unseen nebulae imply values of t\r 
< 10 ~ 5 (jFrail fc Scharringhausenl 119971 iGaensler et al I 
2000). In the case of an undetected, extended radio neb- 
ula surrounding PSR B1046— 58 and using a value of r\ R 
= 2x 10 -4 , the upper limit on the radio flux and required 
surface bright ness imply a large radius of 8 pc for a circu- 
lar nebula ijStappers et al.119 99'). At a distance of 2.7 kpc 
this represents a radius of ~10'. While the X-ray size of 
PWNe is often found to be up to a few times smaller 
than the radio size due to smaller syn chrotron lifetimes 
in the X-rays (e.s.. iHester et aLll2002ft . we find no evi- 
dence for extended PWN structures on arcminute scales. 
For radio sizes < 30", based on the small angular size of 
the X-ray nebula, the required efficiencies are very low 
at n R < 5xl0~ 7 . 

The origin of PWNe is commonly attributed to the in- 
teraction of a highly relativistic pulsar wind with its sur- 
roundings. At a radius r s from the pulsar, representing 
the point of pressure balance at which the wind is con- 
fined and decelerated, we expect P = E / AnQr 2 c. Here, P 
is the surrounding pressure and f2 < 1 is the filling factor 
of the wind. As the head region of the nebula is resolved 
as an extended structure with Chandra, we suggest that 
its outer radius can represent an upper limit on the loca- 
tion of the wind termination shock, while the body repre- 
sents the expected emission downstream. The shock ra- 
dius in this case is r s < 2" = 0.02 6 d?, 7 pc, very simil ar to 
that found for the Vela pulsar (Hclf and et al .11200 1[) and 
an order of magnitude smaller than those of much more 
energ etic pulsars such as the Crab and PSR B15 09— 58 
(e.g.. iWeisskonf et alJl200(H IGaensler et al]l2002ft . The 
required pressure in this case is P > 8.2xl0 -10 d^ 2 ergs 
cm~ 3 (d 2 .7 is the distance to the pulsar in units of 2.7 
kpc). In the region downstream of the shock we expect 
equipartition between the particles and magnetic field to 
be reached, so that B 2 /An = P, where B is the mean 
magnetic field in the nebula. We then estimate a value 
of B > 100 d^i pG. The corresponding synchrotron life- 
time of particles emitting at an energy CkeV ( m units of 

keV) is very small, at i sy nch 40 d^j ^J/' yr. The ve- 
locity that is needed for these particles to reach the edge 
of the nebula within their synchrotron lifetimes is v > 

2,600 d~ 7 /2 eli 2 v km s" 1 > 0.01 d~ 7 /2 e l £ y c. _ 

One problem with the above interpretation is that we 
expect the emission downstream of the shock to be sym- 
metric about the pulsar. A possible explanation for the 
lack of emission to the NW involves Doppler boosting of 
the approaching (in this case SE) component. The ob- 
served flux ratios on either side of the pulsar require an 
intrinsic expansion velocity of vs'm9 > 0.22c, where 6 is 
the inclination of the nebula to the line of sight. Simi- 

9 IGott hclf (2004) has later suggested that this relationship only 
holds for pulsars with _E>4x 10 36 ergs s — 1 and without a bow-shock 
morphology. 



lar velocities are fo und in the structur es of other PWNe, 
such as the Crab (Hes ter et al.ll2002T) . An intrinsically 
asymmetric pulsar wind could also be responsible for the 
observed characteristics. 

The emission in the body of the nebula could also rep- 
resent a collimated outflow, or jet, along the pulsar spin 
axis, as observed in other systems such as the C rab and 
Vela llWeisskopf et al.ll200(il IHelfand et alJl200l . In this 
case, we can estimate the minimum magnetic field and 
energy needed to support the jet against the surrounding 
pressure using equipartiti on arguments (e.s.. iPacholczvld 
119701 iSeward et all 119841) . We assume the jet to be a 
cylinder with length lj = 6.5" = 0.085(^.7 pc and radius 
rj = 1.5" = 0.020(^2. 7 pc, resulting in an emitting volume 
of V ~ 3xl0 51 c?2 7 cm 3 . From the results shown in Ta- 
blcs^a-nd^ we infer an unabsorbed jet luminosity in the 
0.5-10.0 keV range of Lj > 2xl0 31 ergs s" 1 (la limit). 
Equipartition arguments (B oc L 2 / 7 V~ 2 / 7 , E oc B 2 \/) 
then imply a magnetic field in the jet of Bj > 42 (1+p) 2 / 7 
cj)- 2 / 7 pG and an energy ofEj > 4.4xl0 41 (1+m) 4/7 </> 3/7 
ergs. Here, p is the ratio of ion to electron energy and 4> is 
the filling factor in the jet. We note that these represent 
very conservative lower limits on Bj and Ej , as the total 
emission from the jet can have a significant contribution 
at lower energies. Using the above limit on Bj, the corre- 
sponding synchrotron lifetime for particles emitting at 2 
keV is tj < 100 (l+p)~ 3 / 7 3 / 7 yr and implies a velocity 
in the jet of Vj > 1,000 (1+p) 3 / 7 0" 3/7 km s _1 . 

For the Crab and Vela, the nebular magnetic field has 
been estimated to be ~100 's pG (e.g.. lAharonian et alJ 
120041 iPavlov et all l2001a|) . Similar fiel ds have also 
been found for the large-scale jet in Vela l)Pavlov et alJ 
2003|. These are consistent with our derived fields above, 
suggesting the presence of similar confining pressures. 
The size and energetics of the nebula surrounding PSR 
B1046— 58 are similar to that of Vela, as expected if 
these properties are also related to the pulsar's spin- 
down energ y E. In contrast , the n ebula surrounding PSR 
B1509— 58 l|Oaensler et al.l l2002t) shows much smaller 
equipartition fields (~8 pG) and larger size, pointing to 
a smaller confining pressure. 

An alternate argument for an asymmetric nebula in- 
volves the confinement of the pulsar w ind due to the 
ram p ressure of a fast-moving pulsar (e.g.. IGaensler et alJ 
2004J) • In t hi s 

case, the body of the nebula would repre- 
sent the expected X-ray tail aligned with the bow shock 
axis in the direction opposite to the pulsar proper mo- 
tion. The implied position angle in the plane of the sky 
of the proper motion would be ~310° (measured east of 
north). Recent Chandra observations of the Geminga 
pulsar revealed a somewhat similar morpholo gy for its 
PWN l|Pavlov et a,l .1120051 Ide Luca et a,l 1 l2T)06h which is 
likely to arise due the high space velocity of the pulsar. 

In the case of a high space velocity, the confin- 
ing pressure for the nebula is given by P = pV 2 , 
where p is the ambient density and V is the pulsar's 
speed. For cosmic abundances and a number den- 
sity no in the ambient medium we get p = 1. 37nnm.fr 
= 2.3xl0~ 24 no g cm" 3 , where no is in units of 1 
cm" 3 . Our above estimate for P results in a velocity 

°f V ^ 190 n o ^ 2 km s _1 . Such a velocity is consis- 
tent with estimates for other bow-shock nebulae (e.g., 
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iChatteriee fe Corde s 2002) an d the expected distribution 
of pulsar velocities at birth l|Arzoumanian et alJ 119941 
iFaucher-Giguere fe Kaspill2005|) . The proper motion for 
the pulsar in this case would be > 0'.'015 d^\ yr" 1 and 
could be detected by radio interferometry. We note that 
a bow-shock interpretation for the X-ray morphology 
would differ from the static PWN interpretation favoured 
to explain the non-de tection of an extended radio nebula 
ijStappers et alJll999[) . A consistent bow-shock interpre- 
tation in both radio and X-rays is possible if the unde- 
tected radio nebula is unresolved with r s ^ rac iio < 10" and 
has a low efficiency of r/R ~ I0~ 6 , which results in es- 
timates from the ra dio observations of V ~ 250 km s" 1 
and no ~ 0.1 cm -3 (Stap pers et al.lll999D . 

As suggested above, the structures observed in the 
body of the nebula can represent emission due to par- 
ticles left behind after the passage of the pulsar for a 
large enough velocity. For V > 190 km s -1 , the time 
taken by the pulsar to travel the observed distance is 
< 540 d\ 7 yr. Synchrotron lifetimes matching this limit 
require magnetic fields of B > 20 ^G at 1 keV. This is 
consistent with the above equipartition field expected to 
hold in the vicinity of the pulsar given the uncertainties 
in our estimates above. The emission from the body and 
the north clump could also contain contributions from 
high velocity outflows. The north clump at a radius ~ 
3" from the pulsar "ahead" of the bow is not expected 
to arise due to shocked ambient medium. For exam- 
ple, at a velocity V — 190 d^\ km s _1 and density no 
= 1 cm -3 , the medium is expected to have a tempera- 
ture kT ~ 0.07no 1 keV. The luminosity is expected to 
be Lx ~ AripV' ergs s^ 1 ~ 10 28 nQd 3 7 egs s _1 , where 
A ~ I.3xI0~ 24 n o X ^ 2 ergs cm 3 s _1 is the cooling func- 
tion at kT ~ Omn^ 1 keV and V - 7xl0 51 d 3 , 7 cm 3 
is the volume of the shocked gas within 3" of the pul- 
sar (s ee, e.g.. iRvbicki fe Lightmanlfl979l: iGaensler et al.l 
2004) . The implied absorbed flux at Earth is then much 
smaller than the measured values. Additional observa- 
tions can be used to confirm the presence of the north 
clump and its association to the system. 

6.3. PSR B1046-58 and 2EG J1049-5827 

Although not confirmed, an association between 
PSR B1046-58 and 2EG J1049-5827 has been pro- 
posed in light of, among other things, their posi- 
tional coincidence, possible 7-ray pulsations, an d spec- 
tral a nd energetic prop erties of the 7-ray p ulsar l|Fierrol 
IT9951 IPivovaroff etalll2000at iKasni et alJl2000|) . 2EG 
J1049— 5827 was found to have a pulsed 7-ray flux above 
400 MeV of (2.5±0.6)xl0- 10 ergs cm~ 2 s -1 , implying 
a 7-ray efficiency of (9 ±2)xl0~ 3 (for 1 s r beaming and 
a distance of 2.7 kpc. iKaspi et al.l 12000). Unresolved, 
impulsed 7-ray emission at the pulsar position was also 
found. 

An extrapolation of the observed X-ray spectrum for 
PSR Bf 046-58 and its PWN (see Table HJ) to energies 
>400 MeV predicts a flux of 0.7xI0~ 13 ergs cm -2 s _1 , 
three orders of magnitude lower than that derived for the 
EGRET source. Th is is true for other wel l-established 
7-ray pulsars (e.g.. iThompson et al.llT999l) . Observa- 
tions at higher X-ray energies (e.g., deep exposures us- 
ing RXTE, Suzaku) can help determine whether PSR 
B1046— 58 is indeed a 7-ray pulsar by detecting possible 



high X-ray energy pulsations and spectral characteristics 
consistent with those of 2EG J1049-5827. 

We note the presence of a particularly hard X-ray 
source visible in the Chandra data, ^3.'3 from the pulsar 
position at coordinates a 20 oo=10 ?i 48"W74(±0 1 rl5) and 
&.nnn=-58 °28'46"(±2'/4). The sou rce is coincident with 
"Src 4" in IPivovaroff et alJ l|2000aj) : see also Figure El A 
total of 164±15 counts were detected from this source 
and its spectrum is well described by a non-thermal 
power-law model with N H = (0.7±°; 4 ) x f 22 cm~ 2 , V 

= l.lto.2 and unabsorbed flux in the 0.5-10.0 keV of 
(1.2tg;|jxl0 _13 ergs s" 1 cm" 2 ( X 2 = 0.45 for 5 d.o.f., 
la errors). No counterpart for this source was found 
at other wavelengths and it lies within the 99% confi- 
dence error of 2EG J1049-5827. This suggests that the 
observed 7-ray emission could have contributions from 
more than one source. Higher spatial resolution 7-ray 
observations with future missions may help to address 
this possibility. 

We also find that the emission from the nebula dom- 
inates at high X-ray energies (see Table ^) over that of 
the pulsar. This argues for the unpulsed 7-ray emission 
associated with 2EG J1049— 5827 to also arise, at least 
in part, from t he nebula, as prop osed for the Vela and 
Crab pulsars iNolan et alJ Il993t iKanbach etHI 119941 
iFierro et alJll998j) . 

6.4. Non-detection of a Supernova Remnant 

As we do not detect emission from a supernova rem- 
nant ( SNR) in the field cont aining PSR Bf046-58, fol- 
lowing IGaensler et alJ l)2003|) we can estimate the flux 
from a possible unseen remnant by scaling the emission 
observed from the Vela SNR. The emission from this rem- 
nant can be described by a two-component Raymond- 
Smith pl asma with temperatures fcXi ^0.15 keV and kTi 
~ 1 keV l)Lu fc Aschenbachl2 000). The total unabsorbed 
luminosity in the 0.1-2.4 keV range is 2.2 xlO 35 ergs s _1 
and the cooler component has a volume-integrated emis- 
sion measure ~10 times higher than the hot component. 
Scaling to a distance of 2.7 kpc, the unabsorbed flux den- 
sities for each component would be fi — 2.3xl0 -10 ergs 
s _1 cm -2 and fa = 2.2xf0~ n ergs s _1 cm~ 2 , respec- 
tively. 

The expected ACIS-S 10 count rate predicted by the 
WebPIMMS tool 11 is <2.6 counts s" 1 (on-axis) for ab- 
sorbing columns >5xl0 21 cm" 2 . Scaling the size o f the 
Vela SNR (-8° at 300 pc. ILu fc Aschenbachl 1200(1) . we 
would expect a shell with a radius of ~26' at 2.7 kpc. 
This size is larger than the ACIS-S (and EPIC-MOS) 
field of view. However, parts of the remnant might be 
visible in the present data given the uncertainties in 
the spatial distribution of the remnant and direction of 
motion of the pulsar. Assuming emission areas for the 
remnant proportional to having a shell thickness >20% 
of its radius, the expected surface brightness density is 
<0.8xI0~ 6 counts s" 1 arcsec -2 . This is lower than the 
ACIS-S background level during the observation, esti- 
mated to be ~2— 6xl0 -6 counts s _1 arcsec -2 (account- 

10 The EPIC-MOS data provide less constraining limits for this 
emission due to their lower sensitivity and higher background lev- 
els. The EPIC-PN data had a much smaller exposure time and 
limited field of view. 

http:/ /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 
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ing for vignetting at large off- axis angles). 

It is then possible that a remnant with similar proper- 
ties to that of the Vela remnant might be present and it 
is undetected by these observations. The lack of radio or 
X-ray emission from the SNR has also been attributed to 
an initial fa st exp ansion into a low-density cavity (e.g., 
iBraun et alJ IT989). This expansion is followed by a colli- 
sion with a dense surrounding shell of stellar wind mate- 
rial from the progenitor where the emission fades rapidly 
and energy is dissipated through radiative shocks. 

7. CONCLUSIONS 

We have detected an asymmetric, arc-second scale 
PWN surrounding the young, energetic PSR B1046— 58. 
The overall emission from the pulsar and PWN is 
best described by a non-thermal power-law model with 
r=1.7^g'2 and a low X-ray luminosity of ~10 32 ergs s _1 
in the 0.5—10.0 keV range. The brightest emission region 
in the PWN is coincident with the radio coordinates of 
the pulsar. It is resolved as an extended structure with 
Chandra and we suggest that it can represent an upper 
limit to the location of the wind termination shock. The 
additional, asymmetric structures can represent nebular 
emission downstream of the shock or collimated outflow 
from the pulsar. The implied equipartition magnetic 
fields are > 40—100 /iG. The overall size and energet- 
ics of the system in this case are very similar to those of 
the Vela pulsar. Instead, if the observed asymmetry in 
the nebula were to arise due to a large space velocity for 
the pulsar we estimate a value of > 190 n^ 1 ^ 2 km s^ 1 . In 
this case, we would favor the presence of arcsecond-scale 
nebulae with low efficiencies in both X-ray and radio. 
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Our spatially resolved analysis also hints at the presence 
of softer emission from the pulsar than from the rest of 
the nebula. In the case that thermal emission from the 
whole surface of the star is present, we derive an upper 
limit for the temperature of <1.4xl0 6 K, which is not 
constraining on cooling models of neutron stars. In the 
case of a small hot spot on the surface, the implied tem- 
perature is >2.5xl0 6 K and would be consistent with 
those seen for other pulsars. 

The results shown here demonstrate the need for high- 
resolution, high-sensitivity observations in order study 
the wide range of structures associated with rotation- 
powered pulsars. While the origin of these structures 
cannot be unambiguously determined in the case of PSR 
B1046— 58, they can be broadly understood using current 
theories for the production PWNe. Additional data and 
theoretical work are needed to reach a consistent picture 
of this interesting phenomenon in neutron star physics. 
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TABLE 1 

Spatial analysis of PSR B1046-58 and its PWN with Chandra 





0.5-2.0 keV 


2.0-10.0 keV 


HR = 


S/N 


Region 


counts (S) 


counts (H) 


(S-H)/(S+H) 


(0.5-10.0 keV) 


Pulsar (1" radius a ) 


43 ± 8 


24 ± 6 


0.27 ± 0.10 


8.5cr 


PWN, all (- Pulsar) 


46 ± 9 


70 ± 11 


-0.21 ± 0.06 


9.1fT 


PWN, "head" (- Pulsar) 


12 ± 5 


31 ± 7 


-0.44 ± 0.21 


6.1(T 


PWN, "body" 


22 ± 6 


34 ± 7 


-0.22 ± 0.10 


9.1(7 


"North clump" 


5 ± 2 


5 ± 2 


0.00 ± 0.02 


2.9cr 


Pulsar + PWN 


89 ± 11 


95 ± 12 


-0.03 ± 0.01 


12.1(7 



a See JJfor detailed discussion. 



TABLE 2 

Power-law Fit to combined emission 
from PSR B1046-58 and its PWN a 



Parameter Value (±1(7) 



N„ (10 22 cm" 2 ) 


9+ ' 4 
u.a_ 2 


r 


i 7+0-4 
1 -'-0.2 


X 2 (dof) 


79(68) 


fabs b (ergs s" 1 cm -2 ) 


0.7t^xl0- 13 


funabs b (ergs s" 1 cm" 2 ) 


I.O+qjXIO -13 


Lx c (ergs s _1 ) 


0.9l°;6 x i0 32 



a Simultaneous fit to both Chandra and 
XMM- Newton spectra. 

b Absorbed and unabsorbed X-ray fluxes, 
f abs and f un abs, in the 0.5-10.0 keV range. 

c Unabsorbed X-ray luminosity in the 
0.5—10.0 keV range for a distance of 2.7 kpc. 
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